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Avian coronavirus infectious bronchitis virus (IBV) poses economic threat to the poultry industry 
worldwide. Pathogenic IBV 3575/08 was isolated from broilers vaccinated with the attenuated viral 
vaccine derived from a Taiwan strain 2575/98. In this study, extensive investigations were conducted on 
the genome sequences, antigenicity, pathogenicity, and host immune responses of several IBV strains 
in specific-pathogen-free chickens. Sequence analyses revealed that 3575/08 and 2575/98 shared high 
homology in their structural genes, but not in non-structural accessory proteins such as 3a, 3b and 5b. 
Despite a high degree of homology in their spike protein genes, cross neutralization test showed low 
cross protection between 3575/08 and 2575/98, suggesting distinct antigenicity for the two strains. 
Animal challenge experiments exhibited strong respiratory and renal pathogenicity for 3575/08. In 
addition, early and prolonged viral shedding and rapid viral dissemination were observed. Immune 
gene expression profiling by PCR array showed chickens infected with 3575/08 had delayed expression 
of a subset of early innate immune genes, whereas chickens infected with the wild-type or attenuated- 
type 2575/08 revealed quick gene induction and efficient virus control. In summary, this study reveals 
a new IBV strain, which harbors a known local genotype but displays remarkably altered antigenicity, 
pathogenicity and host defenses. 


Avian coronavirus infectious bronchitis virus (IBV) is a member of Coronaviridae family, and the virus is respon¬ 
sible for infectious bronchitis, an acute, highly contagious chicken viral disease. Respiratory signs such as cough¬ 
ing, sneezing, and nasal discharge are the primary symptoms associated with IBV infections, and some of the 
strains also display kidney and oviduct tropisms, causing urogenital symptoms 1-3 . The disease poses economic 
threat to the poultry industry worldwide. More than 50 serotypes of IBV have been documented since the first 
isolation of virus in Massachusetts (Mass), US in the 1930 s 2 . The most commonly detected genotype of IBV in the 
US is Ark-DPI, accounting for 23% to 65% of total isolation per year 4 . While most strains circulating in the US 
cause respiratory symptoms, nephropathogenic QX-like strains are dominant in Europe 2 . In the Far East, nephro- 
pathogenic LX4, LDL as well as proventricular-type QX-like strains have been reported frequently 2 . Presently, 
vaccination is adopted to manage the disease, and Mass type vaccine HI20 is the most widely used vaccine for¬ 
mulation 2,5 . Due to frequent strain variations, however, traditional Mass type vaccines can be ineffective against 
infections by different IBV variants 2,5,6 . Effective disease management thus demands better elucidation of strain 
variations. 

In Taiwan, viral surveillance from 2005 to 2006 in poultry slaughterhouses showed the IBV flock prevalence 
was 17% and the chicken prevalence ranged from 1.8 to 3.7% 7 . Most of the IBV isolates in Taiwan are nephrop¬ 
athogenic strains 8,9 . Also, IB Vs in Taiwan are distinct from other genotypes in the world and can be divided into 
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two subgroups, Taiwan Group I (TW-I) and Taiwan Group II (TW-II) 10,11 . As TW-I strains begin to emerge in 
neighboring areas in recent years 12-14 , there is a rising urgency to thoroughly investigate the TWIBV strains. 

IBV virion has a linear, positive sense and single-stranded genomic RNA. Its 27.6 kb genome contains at least 
ten open reading frames (ORFs), S'-la-lab-S^a-Sb-E-M-Sa-Sb-NW, encoding four structural proteins and sev¬ 
eral non-structural proteins (nsp) 3 . The four structural proteins include spike glycoprotein (S), envelope protein 
(E), membrane glycoprotein (M), and nucleocapsid protein (N). Of the different structural proteins, the S protein, 
which comprises of SI and S2 subunits, is responsible for host cell attachment and virus-cell membrane fusion 
and is the primary antigen target for virus neutralization. Currently, genetic grouping for IBV is based primarily 
on the SI gene 1-3 . Among the genes that encode for non-structural proteins, the ORF la and lab encode two 
polyproteins, which are further cleaved into 15 nsps. These and other non-structural accessory proteins 3a, 3b, 
5a and 5b are associated with RNA replication, virulence, and host-pathogen interactions of IBV 3 . However, the 
functions of these nsps remain largely unknown and warrant further studies. 

Host innate antiviral immunity plays a critical role in viral clearance. In general, toll like receptor (TLR) fam¬ 
ilies, retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated gene 5 (MDA5, encoded by 
IFIH1 gene) are intracellular pattern-recognition receptors (PRRs) that function as viral infection detectors and 
all contribute to type I interferon (IFN) production 15 . As RIG-I is absent in chickens, IFIH1 is a functional com¬ 
pensate 16 . Avian TFRs are different from their mammalian counterparts, lacking TFR9 but possessing TFR1FA, 
poultry-specific TFR15, and a TFR9 homologue TFR21 17 . TFR3 and TFR7 are known to be involved in RNA 
virus recognition, and its stimulation leads to activation of IRF3 and IRF7 (IFN-regulatory factor 3 and 7). These 
factors then signal IFN transcription to trigger IFN-stimulated genes (ISGs) through JAK-STAT pathway to 
express various antiviral proteins 15,17,18 . A global gene expression profiling of chicken kidney has revealed that 
103 immune and inflammatory-related genes are involved in the host defense response during the first week of 
IBV infection 19 . The TFR3/7 and MDA5 signaling pathways and innate immune cytokine were induced after IBV 
infection in chicken tracheas at 3-8 days post IBV M41 infection 20 . In addition, STAT1 expression has also been 
found to be significantly elevated in response to IBV infection in ovo 21 . However, there is a lack of understanding 
of early antiviral immune responses induced by various IBV strains, particularly responses triggered by strains 
with close genetic relationship but distinct pathogenicity 

To control IBV in Taiwan, an embryo-attenuated strain (2575 AT) derived from IBV TW-I strain 2575/98 
wild type (2575 WT) has been developed as a vaccine 8 . These two strains were previously investigated to have 
close antigenicity 8 and limited sequence changes 22 . However, a novel TW-I strain 3575/08 emerged from a broiler 
chicken farm where the 2575 AT vaccine was applied. IBV 3575/08 has been observed to cause high morbidity 
and mortality in affected broilers. In the present study, we compared the viral genome sequences of IBV 3575/08 
and 2575/98 and investigated the infectivity and host immunity induced by 3575/08, 2575 WT, and 2575 AT 
infections. To this end, sequence analysis, neutralization test, experimental infection, and immune-related gene 
profiling were performed. The results show that IBV 3575/08, which possesses heterologous ORF 3a, 3b, and 5b 
sequences to 2575/98, caused a more deadly infection with early and prolonged viral shedding, stronger clinical 
and pathological manifestation, and delayed innate immune responses as compared to 2575 WT and 2575 AT. 

Results 

3575/08 and 2575 WT share sequence homology except in 3a # 3b # and 5b genes. The complete 
genome of IBV 3575/08 was sequenced in this study, and the sequences were submitted to GenBank with the 
accession number KX266757. Excluding the 5' and 3' untranslated regions the entire genome size is 26,704 
nucleotides. The genome comprises ten ORFs in the order of 5 / -la-lab-S-3a-3b-E-M-5a-5b-N-3 / . No significant 
protein truncation was observed. To evaluate the nucleotide sequence identities among 3575/08 and other refer¬ 
ence IBV strains, sequence alignment and phylogenetic analyses on each gene segment were performed (Fig. 1 
and Supplementary Fig. SI). Comparisons of the available full-length genome showed that the sequence identi¬ 
ties between 3575/08 and other IBV strains ranged from 85.0% (with strain BJ) to 94.9% (with strain 2575/98). 
Overall, IBV 3575/08 showed close genetic grouping with 2575/98 and other Taiwan local strains in la, lb, S, E, 
M, 5a and N genes. The sequence similarities of these protein genes between 3575/08 and 2575/98 were 95.6% 
(la), 94.2% (lb), 95.6% (S), 93.3% (E), 92.4% (M), 93.4% (5a), and 96.3% (N) for the respective genes. Further 
comparison of the two spike protein subunits between 3575/08 and 2575/98 showed sequence similarities of 
97.1% and 94.4% in the SI and S2 subunits, respectively. Based on the S gene, 3575/08 was clearly defined as a 
TW-I strain. The sequence identities of S gene between 3575/08 and other analyzed strains ranged from 67.6% 
(with strain Georgia 1998 pass8) to 96.4% (with strain 3468/07). While 3575/08 was closely related to 2575/98 
in four structural protein genes, their homology in ORF 3a (86.1%), 3b (82.1%), and 5b (95.5%) were relatively 
low (4.5% difference is considered significant for the 5b protein since its variation among IBV strains is small). 
According to the phylogenic tree (Fig. 1), IBV 3575/08 and 2575 WT are grouped in different clusters with respect 
to their 3a, 3b and 5b genes. 

Identification of IBV 3575/08 as a unique serotype. To investigate the antigenic serotype of the newly 
emerged TW-I strain 3575/08, cross (two-way) neutralization tests against three major serotypes in Taiwan— 
TW-I, TW-II and Mass type— were performed in specific-pathogen-free (SPF) chicken embryos. The neu¬ 
tralizing titers and the r-values were then analyzed (Table 1 and Supplementary Table SI). With the r-value of 
homologous virus-serum set as 100, the r-value directly correlates with the relatedness of the strains. The results 
showed that all the r-values obtained in this study were below 30, indicating all the tested strains belonged to dif¬ 
ferent serotypes. Particularly, TW-I 3575/08 conferred poor cross neutralization not only to heterogeneous Mass 
type H120 (r-value = 1.1) or TW-II 2296/95 (r-value = 3.4) but also to homologous TW-I 2575/98 (r-value = 3.3), 
supporting IBV 3575/08 as a unique serotype. 
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Figure 1. Phylogenetic analyses of la, lb, S, 3a, 3b, 5a and 5b protein genes. The phylogenetic trees were 
constructed using the MEGA software version 4 by the neighbor-joining method (bootstrapping for 1,000 
replicates with its value >70%). The 3575/08 and 2575/98 strains analyzed in this study are indicated with red 
triangles. (AU: Australia, CN: China, NL: Netherlands, TW: Taiwan, UK: United Kingdom, US: United States, 
WA: Western Africa) 


IBV 3575/08 showed enhanced pathogenicity as compared to 2575 WT and 2575 AT in chickens. 

To investigate the pathogenicity of IBV 3575/08, 2575/98 WT and 2575 AT, SPF chickens were intranasally 
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IBV strain 

3575/08 (TW-I) 

2575/98 (TW-I) 

2296/95 (TW-II) 

H120 (Mass) 

3575/08 (TW-I) 

100 




2575/98 (TW-I) 

3.3 

100 



2296/95 (TW-II) 

3.4 

12 

100 


HI20 (Mass) 

1.1 

2.6 

12.8 

100 


Table 1. R-value of cross neutralization test with TW-I strains (3575/08, 2575/98), TW-II strain (2296/95) 
and Mass type (HI20). 


challenged with 10 6 EID 50 of the three strains and monitored for survival and clinical signs on a daily basis. In 
1-day-old chicks during a 21-day observation period, IBV 3575/08 infection resulted in the lowest survival rate 
(16.7%, 1/6) as compared to 2575 WT (66.7%, 4/6), 2575 AT (100%, 6/6), and PBS control (100%, 6/6) (Fig. 2A). 
The earliest death from 3575/08 or 2575 WT-infected group was found on day 4 post infection (dpi). 

In addition, clinicopathological assessment for the three IBV strains was performed. Following intranasal 
challenge in 2-week-old chickens with 10 6 EID 50 of the different IBV strains, the clinical scores (Fig. 2B,C) were 
recorded and analyzed (n = 5 per group, see Methods for detailed clinical score interpretation). It was observed 
that beginning at 2 dpi, chickens infected with 2575 WT exhibited moderate tracheal rales (score =1), persisting 
through 3 dpi (score = 1) and resolved at 4 dpi (score = 0.2) (Fig. 2B). IBV 3575/08-infected chickens showed 
similar but more severe clinical signs at 2 dpi (score = 1.4) and 3 dpi (score = 1.2). At 4 dpi, more than half of 
3575/08-infected chickens still suffered from mild respiratory symptoms (score = 0.6), and clinical signs receded 
after 5 dpi. In contrast, IBV 2575 AT-infected and PBS control groups did not show any clinical manifestations 
throughout the observation period. Collectively, over a 7-day observation period, the averaged clinical scores of 
IBV 3575/08-, 2575 WT-, 2575 AT-infected chickens were 0.49, 0.32, and 0, respectively (Fig. 2C). Overall, IBV 
3575/08 infection caused stronger clinical manifestations in chickens than 2575 WT and 2575 AT. 

The pathological lesions of experimentally infected chickens were further evaluated at 7 dpi, and representa¬ 
tive pictures are shown in Fig. 3. In the trachea epithelium, normal goblet cells and ciliated cells were seen in the 
PBS control group (Fig. 3D). However, the IBV 3575/08 infected chickens presented pathological lesions such 
as desquamation of ciliated cells, infiltration of heterophils, and lymphocytes and epithelial hyperplasia in the 
trachea (Fig. 3A, arrowed). Fesions in the 2575 WT (Fig. 3B) and 2575 AT (Fig. 3C) groups were similar but less 
severe as compared to the 3575/08 group. Upon observing the kidneys, the control group showed normal renal 
tubular epithelial cell with homogeneous plasma and distinctive renal tubular cavities (Fig. 3H). In the kidneys 
of IBV 3575/08-infected chickens, however, renal tubule dilation, necrotic tubular cells, and urate in the tubular 
cavity were observed (Fig. 3E, arrowed). Vesicular degeneration was also observed in the renal tubular epithelial 
cells from the 2575 WT group (Fig. 3F, arrowed). On the other hand, the 2575 AT group showed normal kidney 
patterns similar to the control group (Fig. 3G). The pathological findings indicate that both IBV 3575/08 and 
2575 WT are respiratory and nephropathogenic strains. IBV 3575/08 showed higher virulence as it inflicted more 
severe trachea and renal tissue damages in chickens than the 2575 WT strain. No virulence in the kidney was 
observed for IBV 2575 AT. 

More rapid and prolonged viral dissemination of IBV 3575/08 in chickens. To further investigate 
the viral dissemination of IBV 3575/08, 2575 WT, and 2575 AT in chickens, swabs and tissues were harvested 
from virus-infected chickens (10 6 EID 50 , intranasal) at various time points. Viral detection was carried out by 
reverse transcription polymerase chain reaction (RT-PCR). As indicated in Table 2, 100% of the throat swabs 
from all groups were detected positive at 1, 2, 3, 7, 14 dpi. At 21 dpi, both the 2575 WT and 2575 AT groups 
showed a decreased positive rate at 66.7% (2/3), whereas a 100% (3/3) positive rate remained for the 3575/08 
group. More importantly, for the cloacal swabs, virus was detected from 3575/08-infected chickens as early as 
1 dpi (60%, 3/5), whereas the earliest time points for 2575 WT and 2575 AT were 2 dpi (40%, 2/5) and 3 dpi (20%, 
1/5), respectively. There was no detectable level of viruses in any of the swab samples from the PBS control. 

To examine viral distributions, the lungs and kidneys, two major organs associated with IBV tropism, were 
harvested at 3, 7, 21 dpi. As indicated in Table 3, whereas all of the infected chickens from the 3575/08 or 2575 
WT groups were detected positive in the lungs at 3 and 7 dpi, only 2/3 and 1/2 chickens were detected positive in 
the 2575 AT group. As for the kidneys, virus was detected from the 3575/08-infected chickens as early as 3 dpi, 
whereas no virus was detected from the 2575 WT group until 7 dpi. In addition, IBV 3575/08 was detected in the 
trachea, proventriculus, bursa, and oviduct from experimentally infected chickens at 7 dpi (data not shown). In 
the 2575 AT group, no animal was detected positive in the kidneys across the indicated time points. By 21 dpi, 
virus was cleared from the lungs and kidneys in every group. These results show that IBV 3575/08 causes more 
rapid and prolonged dissemination and shedding in chickens as compared to the other two strains. The results 
also highlight the nephropathogenic nature of 3575/08 and 2575 WT in chickens. 

3575/08 infection elicited a similar level of antibody titer as 2575 WT in chickens. To examine 
the antibody response elicited by IBV infection in chickens, sera were collected from the 3575/08-, 2575 WT-, and 
2575 AT-infected groups at 14 and 21 dpi. Antibody titers were analyzed by an IBV antigen coated-EFISA, and 
3575/08 and 2575 WT infections elicited similar mean titer levels of 6,547 and 6,870 at 21 dpi respectively (Fig. 4). 
The antibody titer elicited by 2575 AT infection was lower (3,294) at 21 dpi. Similar titer results were obtained 
from the different infection samples at 14 dpi. 


SCIENTIFIC REPORTS | 6:37725 | DOI: 10.1038/srep37725 


4 














www.nature.com/scientificreports/ 




3575/08 
2575 WT 
2575 AT 
Control 


Days post infection 


C 

0.8- 
g 0.6- 

o 

</> 

ro 0.4- 
o 

c 

0 0.2- 


ns _ **_ 

O 

oo 

A 

oo sags 


0 . 0- 1 - 1 - 

3575/08 


—,-ehhb- 

2575 WT 2575 AT 
IBV strain 


Figure 2. Percentage survival and clinical score of IBV infected chickens. (A) Survival rate of 1-day-old 
chicks after intranasal infection with IBV 3575/08, 2575 WT, 2575 AT. Chickens receiving PBS were served 
as control, (n = 6 per group) (B) Following intranasal infections in 2-week-old chickens, the daily averages of 
observed clinical scores within 6 days were recorded. Data was presented as means d= s.e.m. (n = 5 per group) 
(C) Mean clinical scores of each group during the observation period are indicated (n = 5 per group). See 
Methods for detailed clinical score interpretation, ns: non-significant. **P< 0.01. ***P< 0.001. 


Antiviral immune response was delayed in IBV 3575/08-infected chickens. Next, to dissect the 
antiviral gene induction in chickens in response to IBV 3575/08,2575 WT and 2575 AT infections, we examined 
a panel of immune response gene expression in the lung at 24 and 48 hours after the IBV infections using a quan¬ 
titative PCR array kit. The fold changes in the expression of 84 genes from infected chickens were normalized to 
that of naive chickens (Fig. 5A, also see Supplementary Table S3 for statistical data and Supplementary Fig. S2 
for the Venn diagram). At 24 hours post infection (hpi), 20 (23.8%) and 33 (39.3%) of the genes were induced by 
more than 3-fold over the naive control in the 2575 WT- and 2575 AT-infected chickens respectively, whereas 
only 5 (6%) of the 84 genes were upregulated in the 3575/08-infected chickens. At 48 hpi, 16 (19%), 3 (3.6%), 
and 3 (3.6%) genes were induced by more than 3-fold over the naive group in the 3575/08-, 2575 WT-, and 2575 
AT-infected chicken, respectively Overall, a robust immune-related gene expression was rapidly induced in the 
2575 WT- and 2575 AT-infected chickens at 24hpi, and a 90% overlap was observed among the upregulated genes 
between the two groups. These results identified the critical components in early host defense against IBV 2575 
WT and 2575 AT. At the same time point, 3575/08-infected chickens showed heightened expression of only 5 
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Figure 3. Pathological examination of trachea and kidney from IBV infected chickens. The pathological 
lesions of experimentally infected chickens were evaluated at 7 dpi, and representative pictures were shown. 
Panels D and H refer to control tissues. IBV 3575/08 infected chickens presented pathological lesions such as 
desquamation of ciliated cells, infiltration of heterophils, and lymphocytes and epithelial hyperplasia in the 
trachea (A, arrowed). The lesions in 2575 WT (B) and 2575 AT (C) groups were similar to the 3575/08 group 
but milder. For IBV 3575/08-infected chicken kidney, renal tubule dilation, necrotic tubular cells, and urate in 
the tubular cavity were observed (E, arrowed). Vesicular degeneration was found in renal tubular epithelial cells 
from the 2575 WT group (F, arrowed). For the 2575 AT group, the renal cells were normal as the control group 
(G). 


innate immunity genes (TLR1LA, TLR3, IRF7, STAT1, and C3) (Fig. 5B). 6 other genes upregulated by 2575 WT 
and AT at 24hpi, including TLR4, MyD88, IFNA3, IFNGR1, CD28, and CATH2, showed delayed stimulation fol¬ 
lowing the 3575/08 infection (Fig. 5C). In addition, 10 other genes that were not enhanced in either the 2575 WT 
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IBV strain 

Throat swab % IBV positive 

Id 

2d 

3d 

7d 

14d 

21 d 

3575/08 

100% (5/5) 

100% (5/5) 

100% (5/5) 

100% (5/5) 

100% (3/3) 

100% (3/3) 

2575 WT 

100% (5/5) 

100% (5/5) 

100% (5/5) 

100% (5/5) 

100% (3/3) 

66.7% (2/3) 

2575 AT 

100% (5/5) 

100% (5/5) 

100% (5/5) 

100% (5/5) 

100% (3/3) 

66.7% (2/3) 

Cloacal swab % IBV positive 

3575/08 

60% (3/5) 

60% (3/5) 

100% (5/5) 

100% (5/5) 

100% (3/3) 

100% (3/3) 

2575 WT 

0% (0/5) 

40% (2/5) 

100% (5/5) 

100% (5/5) 

100% (3/3) 

100% (3/3) 

2575 AT 

0% (0/5) 

0% (0/5) 

20% (1/5) 

80% (4/5) 

100% (3/3) 

100% (3/3) 


Table 2. Rate of IBV detection by RT-PCR in chicken throat and cloacal swabs after infections. 


IBV strain 

Lung 

Kidney 

3d 

7d 

21 d 

3d 

7 d 

21 d 

3575/08 

3/3 

2/2 

0/3 

1/3 

2/2 

0/3 

2575 WT 

3/3 

2/2 

0/3 

0/3 

2/2 

0/3 

2575 AT 

2/3 

1/2 

0/3 

0/3 

0/2 

0/3 


Table 3. IBV detection by RT-PCR in chicken lungs and kidneys at 3,7, and 21 day post infection. 


15000n 


■ 14 dpi 
□ 21 dpi 



3575/08 2575 WT 2575 AT 


IBV strain 


Figure 4. Serum antibodies of IBV infected chickens. Chicken serum was collected at 14 and 21 dpi, 
respectively. ELISA was performed to evaluate the antibody level against IBV. Bars represent means d= s.e.m. 
(n = 5 per group). 


or the 2575 AT group were found to be upregulated in the 3575/08 group at 48hpi. The finding suggests that IBV 
3575/08 may have evaded the initial virus sensing and delayed the signaling pathways for early antiviral defense. 
In 2575 WT-infected chickens, the genes upregulated at 48hpi (TLR1LA, TLR3, and CATH2) overlapped with 
those at 24hpi. Similarly, the genes upregulated by 2575 AT at 48hpi (TLR1LA, IFNA3, and CD28) also over¬ 
lapped with those observed at the earlier time point. In addition, quantification of TLR3 and TLR4 expression 
levels by regular real-time RT-PCR was carried out to validate the results from the PCR array The expression 
patterns of these two genes were consistent between the real-time PCR and the PCR array analysis, showing acti¬ 
vation of TLR3 but not TLR4 level at 24hpi in the 3575/08 group (Fig. 5D). 

Discussion 

The spread of TW type IBV in Taiwan has been attributed to poor protection by commercial Mass type vaccines 23 . 
An effective TW type vaccine strain 2575 AT has been previously developed through serial passages in SPF 
embryonated eggs to control IBV infection in Taiwan. However, a new IBV strain 3575/08 emerged from 2575 
AT-vaccinated broiler chickens and caused high morbidity and mortality To characterize the new strain 3575/08, 
comparative analyses on viral genotype, serotype and pathogenicity were conducted. 

Pairwise sequence comparison and phylogenetic analysis revealed that 3575/08 and 2575/98 share high 
homology in most of the structural protein genes, including the SI gene that are commonly applied for IBV strain 
genotyping. Therefore, IBV 3575/08 was grouped in the TW-I genotype with 2575/98 and other strains. It has 
been reported that the serotype of IBV is associated with the hypervariable region (HVR) 1 in the N-terminus 
(residues 45-114) of the SI gene 24 , which overlaps with one of the antigenic sites located in the spike protein of 
IBV 25 . Despite high sequence similarities between 3575/08 and 2575/98 were found in the spike protein gene 
(95.6%), SI subunit gene (97.1%) and HVR 1 (94.8%), a distinct serotype of 3575/08 was clearly established in 
this study through a cross neutralization test. Similar findings were previously reported in UK showing that IBV 
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Figure 5. Comparison of immune related genes in IBV infected chickens. (A) Profiles of immune related 
gene expression in the lungs after infection by different strains as evaluated by RT 2 Profiler Chicken Innate & 
Adaptive Immune Responses PCR Array at 24 and 48 hpi. Gene expression of transcripts with a difference of 
more than 3-fold over naive was illustrated (black, upregulated; gray, no change; white, downregulated). 

(B) Gene induction in the lungs of IBV 3575/08,2575 WT, and 2575 AT-infected chickens at 24 hpi was 
expressed as fold-increase over naive. (C) Gene induction in the lungs of IBV 3575/08 and 2575 WT-infected 
chickens was expressed as fold-increase over naive at 24 or 48 hpi. (D) Expression level of TLR3 and TLR4 
at 24 hpi in was validated by real-time RT-PCR. Fold changes were compared by ANOVA followed by 
Dunnett multiple comparisons tests. Dotted lines denote a 3-fold change over naive chickens. Bars represent 
means is.e.m. (n = 3 per group). *P< 0.05, **P< 0.01. 
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strains of three distinct serotypes possess as little as 2% of amino acid variations in region 19-122 and 251-347 
of the SI subunit 26 . The variations led to failure of cross-protection 27 , supporting that small sequence differences 
can result in a serotype change. In addition, it has been reported that small sequence changes in the S2 subunit 
may alter interactions with the SI subunit, which could in turn affect the conformation of the SI subunit and 
serotype-specific epitopes 28 . A previous study on murine hepatitis virus (MHV), a group II coronavirus, also 
showed that a single amino acid substitution in region 1109-1116 of S2 subunit can confer resistance to neu¬ 
tralization by SI subunit-specific monoclonal antibody 29 . Therefore, the gene variations in the S2 subunit (5.6%) 
between 3575/08 and 2575/98 may have a significant contribution to the serotype difference. Future studies are 
warranted to further examine the determinants of antigenicity. 

In addition to the new serotype, we also found that IBV 3575/08 exhibited higher mortality and more 
profound clinical signs as compared to 2575 WT and 2575 AT through experimental challenges in chickens. 
Following intranasal infection, IBV 3575/08 had a very short latent period. The viral shedding from cloaca was 
detected as early as 24hpi and lasted for 21 days, and a rapid viral dissemination to renal tissues (3 dpi) was 
observed. As it has been documented that spike proteins play a major role in determining viruses’ cell tropism 
and pathogenicity 30 , the increased IBV 3575/08 virulence is a notable finding given the high degree of homology 
between the spike protein genes between IBV 3575/08 and 2575 WT. 

A previous study has shown that different IBV strains can cause differing IFN levels at various organs and 
time points in chickens 31 . We hypothesized that a different innate immune response may account for the altered 
pathogenicity of IBV 3575/08. Therefore, the early immune-related gene expression in the lung following 3575/08, 
2575 WT and 2575 AT infections in chickens was analyzed by PCR arrays. In contrast to the 20 and 33 genes 
activated in the 2575 WT- and 2575 AT-infected chickens, only a small subset of the genes was triggered in the 
3575/08-infected chickens at 24hpi, and a delayed activation was observed at 48hpi. This result suggests that IBV 
3575/08 has the potential to evade the initial virus sensing and suppress the first wave of host defense, leading 
to enhanced pathogenicity observed in this study. Many coronaviruses have the ability to evade detection by the 
host’s immune system and disrupt the host immune response through multiple mechanisms 32-34 . It has been 
reported that small accessary proteins 3a and 3b of IBV can delay the onset of IFN response by interfering with 
the PRR sensing without disrupting the signaling pathways 35 . Previous studies have also demonstrated that 3a 
and 3b accessory proteins of IBV are not essential for replication but can modulate the IFN-(3 response at the 
transcriptional and translational level 36,37 . In addition, 5b accessory protein of IBV has been found to shut off 
host immune response and limit interferon production 38 . It has been shown that a natural recombinant strain 
containing partial nsp regions from a low pathogenic HI20 strain in the frame of a nephropathogenic strain 
possesses attenuated pathogenicity, which suggests strong associations between nsps and IBV virulence 39 . Studies 
on other coronaviruses such as MHV and severe acute respiratory syndrome coronavirus (SARS-CoV) have also 
shown that nsps are associated with viral replication and evasion of host innate immune response via viral RNA 
modification 40,41 , host mRNA degradation, and IFN antagonist competition 42 . Therefore, given the low sequence 
homology in the accessory proteins 3a, 3b, and 5b between IBV 3575/08 and 2575 WT, we conjecture that these 
accessary proteins contribute to the blockade of PRR expressions, delayed immune response, and enhanced 
virulence of IBV 3575/08 in the present study. Future studies are warranted to further dissect the roles of IBV 
3575/08’s nsps and accessary proteins in modulating host’s antiviral defenses. An in vitro infection model with 
detailed functional analysis of nsps would shed clearer insights on the host-virus interactions. 

Also of note is that the number of PRR genes activated at the early time point correlated with IBV virulence in 
the study. For instance, IBV 2575 AT infection activated multiple PRR genes, including IFIH1, TLR1LA, TLR2-2, 
TLR3, TLR4, TLR7, and TLR15. In addition to the two primary sensors for IBV, TLR3 and TLR7 43 , other PRR 
molecules have been studied for their roles in coronavirus infections. In previous studies, SARS-CoV and MHV 
infections have been shown to trigger IFIH1 expression for MDA5 signaling 44,45 . It has also been reported that 
TLR4 expression can be induced by SARS-CoV infection 46 . The robust expression of PRR molecules triggered by 
IBV 2575 AT sheds light on the molecular mechanism behind the strain’s attenuated virulence and pathogenicity. 
In contrast, the lower number of stimulated PRR genes triggered by IBV 3575/08 and 2575 WT serve to explain 
these viruses’ increased virulence. 

The present study may also offer insight on IBV disease management and vaccine design. Firstly, identifi¬ 
cation of a new pathogenic IBV strain with a unique serotype and strong host tropism suggests a rapid viral 
mutation rate in chicken hosts. This observation raises the need for regular vaccine update to keep up with the 
mutating viral strains. As the antigenicity or pathogenicity shift may result from genetic changes outside of struc¬ 
tural protein genes, antigen derived from a single epitope (eg. spike protein) may not be sufficient in eliciting 
cross-protection against different serotypes. Multi-epitope vaccine formulations should be considered for better 
disease management. Secondly, the present study elucidates several key innate immune genes and TLR pathways 
(eg. TLR3, TLR4, MYD88, IRF7, IFNA3, etc.) responsible for early restriction of IBV infection. When designing 
an inactivated or subunit vaccine, immunomodulatory adjuvants may be tailored to promote an infection-like 
immune stimulation as observed in the study. Triggering innate immune gene expressions and TLR pathways in 
a virus-like manner may facilitate a robust adaptive immune response against IBV infection. 

In summary, this study compares the sequence, antigenicity, pathogenicity, and immune host defenses among 
IBV strains. The strain 3575/08, which possesses heterologous ORF 3a, 3b, and 5b sequences, exhibits a distinct 
serotype and strong clinicopathological manifestation. This virus can also evade the initial virus sensing and delay 
several early innate immune genes, collectively leading to enhanced pathogenicity in a chicken host. 

Methods 

Eth ics Statement. Experiments involving chickens were conducted at National Taiwan University, under 
an approved Institutional Animal Care and Use Committee (IACUC) protocol (no. NTU-103-EL-3). All animal 
experiments were carried out in accordance with the approved guidelines. 
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Virus. IBV strain 2575/98 was isolated in Changhua, Taiwan in 1998 s , and designated 2575 wild-type (2575 
WT) throughout this study Embryo passage 74 of 2575/98 was obtained as previously described 8 , and designated 
2575 attenuated-type (2575 AT) throughout this study IBV strain 2296/95 was isolated in Taoyuan, Taiwan in 
1995 s . Mass type strain H120 was purchased from Abie Biological Lab (Beit Shemesh, Israel). IBV strain 3575/08 
was isolated from a broiler farm in Yunlin, Taiwan in October 2008. In that affected farm, death was first observed 
in 25-day-old broilers with typical IB disease symptoms. The morbidity was higher than 15%, and the mortal¬ 
ity was approximately 7.5%. All IBVs used in this study were propagated in 10-day-old SPF embryonated eggs 
(Animal Health Research Institute, Tamsui, Taiwan) via an allantoic route as previously described 47 . All the used 
viruses were ensured free of other avian respiratory pathogens, ie. avian influenza virus and Newcastle diseases 
virus contamination by RT-PCR 48 . IBV viral titer (egg infective doses 50%, EID 50 ) was evaluated in 10-day-old 
SPF eggs and calculated by Reed-Muench method 49 . 

Viral genome sequencing and sequence analysis. Viral RNA was extracted from virus-infected allan¬ 
toic fluid using a viral nucleic acid extraction kit (Geneaid Biotech Ltd., Taipei, Taiwan) following the manufactur¬ 
er’s protocol. Reverse transcription was carried out with M-MLV reverse transcriptase (Invitrogen, Carlsbad, CA), 
and genome segments were amplified by PCR using primers listed in Supplementary Table S3. The PCR reaction 
mixture contains 0.5 pi of Phire hot start II DNA polymerase (Thermo, Carlsbad, CA), 5 pi of 5X buffer, 0.5 pi 
of 50 pM forward primer and reverse primer (Tri-I Biotech, Taipei, Taiwan), 3 pi of 2.5 mM dNTPs (GeneTeks 
BioScience, Taipei, Taiwan), 5 pi of cDNA and 10.5 pi of nuclease-free water. For PCR, the initial denaturation was 
at 95 °C for 3 min, followed by 35 cycles of denaturation at 95 °C for 1 min, annealing at 55 °C for 30 s, and exten¬ 
sion at 72 °C for 2 min. The final extension was at 72 °C for 10 min. All the amplified PCR products were analyzed 
in 1% agarose gels electrophoresis. DNA sequencing was conducted by a commercial service (Tri-I Biotech). Each 
nucleotide was determined from at least three identical results generated from separate PCR products. Nucleotide 
sequences were assembled using DNAStar software (DNAStar, Madison, WI). Sequences of the reference IBV 
strains were retrieved from GenBank with the accession numbers listed in Supplementary Table S4. Multiple 
sequence alignments and sequence identity were analyzed with the Clustal W method available in BioEdit 7.2.5 
software. Phylogenetic trees were constructed with the neighbor-joining method using MEGA software version 
4 50 , and the bootstrap values were determined from 1,000 replicates of the original data. 

Antisera production and cross neutralization test. Antisera against 3575/08, 2575/98 and 2296/95 
were prepared in 3-week-old SPF chickens as previously described 24 . The antiserum against Mass-type H120 was 
purchased from Charles River Laboratories (North Franklin, CT). Cross neutralization test was performed with 
the protocol previously described 24 . Briefly, four-fold diluted sera were incubated with the same volume of 100 
EID 50 homologous or heterologous virus at room temperature for an hour, respectively The virus-serum mixtures 
were then inoculated into the allantoic cavity of 10-day-old SPF embryonated eggs. Seven days after inoculation, 
the eggs were opened and examined for typical lesions caused by IBV infection (embryo dwarfing or malforma¬ 
tion). The neutralizing titer of each serum against homologous or heterologous virus was determined by the last 
serum dilution that protected 50% of the embryo and calculated by Reed-Muench method 49 . In addition, the 
cross neutralization r-values between strains were calculated by the method described by Archetti and Horsfall 51 . 
Antigenic (serotype) difference between two given strains was defined as follows, when r: 70-100%, same sero¬ 
type; r: 33-70%, different subtype (minor); r: 11-32%, different subtype (major); r; 0-10%, different serotype 52 . 

IBV infections in SPF chickens. For survival evaluation after IBV infection, 1-day-old SPF chicks were 
randomly divided into three groups and intranasally challenged with 10 6 EID 50 of IBV 3575/08, 2575 WT and 
2575 AT, respectively (n = 6 per group). Another six chickens received PBS as naive controls. For the clinico- 
pathological assessment, 2-week-old SPF chickens were used. Following intranasal challenge with 10 6 EID 50 of 
IBV 3575/08, 2575 WT and 2575 AT, chickens were monitored on a daily basis. The clinical scores of IBV were 
interpreted according to the methods described by Avellaneda et al 53 . The clinical signs were evaluated as: 0 = no 
clinical signs; 1 = slight lacrimation, slight shaking, watering feces or tracheal rales; 2 = lacrimation, presence of 
nasal exudate, depression, water feces, apparent sneezing or cough; 3 = same as 2 but stronger with severe watery 
feces; 4 = death. Throat and cloacal swabs were collected at 1, 2, 3, 7, 14, 21 dpi for virus detection. In addition, 
lungs and kidneys were also harvested from euthanized chickens at 3, 7,21 dpi for virus detection and histopatho- 
logical examination. Sera were collected at 14 and 21 dpi to test the IBV antibody titer. For PCR array analysis, 
three chickens at 24 and 48hpi from each challenged group were euthanized to harvest lungs. Four chickens 
received PBS as naive controls. For all chicken experiments, each virus-challenged group of chickens was reared 
in separated rooms, and utensils were separately handled to avoid any cross contamination. 

Viral detection by RT-PCR. The viral nucleic acid extraction kit (Geneaid) was used to extract the viral RNA 
from swabs samples and tissue homogenates. One step RT-PCR was performed in a total reaction volume of 25 pi 
containing 0.5 pi Phire hot start II DNA polymerase (Thermo), 5 pi 5X buffer, M-MLV RT (200 U/pl) (Invitrogen), 
1 pi RNaseOUT (40 U/pl) (Invitrogen), 0.5 pi 50 pM forward primer (5 / -GGTAGYGGYGTTCCTGATAA-3 / ) and 
reverse primer (5 / -TCATCTTGTCRTCACCAAAA-3 / ) (Tri-I Biotech) for a 618 bp fragment of IBV N gene 
detection, 3 pL 2.5 mM dNTPs (GeneTeks BioScience), 2.5 pi RNA, and 10.5 pi nuclease-free water. RT-PCR was 
performed at 40 °C for 30 min, 98 °C for 3 min, followed by 35 cycles of denaturation at 98 °C for 30 s, annealing 
at 58 °C for 30 s, and extension at 72 °C for 30 s, and the final extension was performed at 72 °C for 3 min. PCR 
products were then analyzed on 2% agarose gels. 

ELISA. A flat-bottomed microplate (Nunc, Denmark) was coated overnight with purified IBV antigen (100 ng/well) 
with coating buffer (15 mM Na 2 C0 3 and 35 mM NaHC0 3 , pH 9.6). The wells were washed three times with PBST 
(0.1% (v/v) Tween 80 in PBS) and blocked with 5% (w/v) skim milk (BD Difco, Sparks, MD) at room temperature 


SCIENTIFIC REPORTS | 6:37725 | DOI: 10.1038/srep37725 


10 




www.nature.com/scientificreports/ 


for 1 h. After blocking, serial diluted chicken serum was added to each well (100 pi/well) and incubated at room 
temperature for 1 h. Following three washes, 100 (Jil of 1:2,000 diluted peroxidase-conjugated affinipure rabbit 
anti-chicken IgG (H + L) (Jackson ImmunoResearch, West Grove, PA) was dispensed into each well and incu¬ 
bated at room temperature for 1 h. After three additional washes, the wells were then incubated with 100 pi of 
SureBlue Reserve TMB Microwell Peroxidase Substrate (KPL, Gaithersburg, MD) for 10 min in the dark. The 
reaction was stopped by the addition of 100 pi of TMB stop solution (KPL). The optical density (OD) at 450 nm 
was read using an automated plate reader (Thermo). 

Histopathological examination. Chicken tracheas, lungs, and kidneys were collected and fixed in 10% 
formalin (Sigma, MO) at room temperature for at least 48 h. Fixed tissues were embedded in paraffin wax, and 
cut into 6 pm sections. Sections were then routinely processed and stained with hematoxylin and eosin (H&E). 
Microscopic lesions resulted from IBV infection were examined by light microscopy 

PCR array. To analyze the innate and adaptive immune responses in chicken after infection, RT 2 Profiler 
Chicken Innate & Adaptive Immune Responses PCR Array (PAGG-052ZG; SABiosciences) was used following 
the manufacturer’s protocol. Total RNA in chicken lung was extracted using TriSolution Reagent Plus (GMbiolab, 
Taipei, Taiwan) following the recommended protocol for tissue samples and quantitated using spectrometry To 
avoid the contaminations of genomic DNA, 1 pg total RNA was treated with 1U DNase I (Thermo) for 30 min at 
37 °C, then 1 pi EDTA was added and incubated 10 min at 60 °C to inactivate the DNase. cDNA was synthesized 
using the RT 2 First Strand kit (SABiosciences), mixed with 2 x RT 2 qPCR Master Mix (SABiosciences), and then 
loaded onto the 384-well PCR array plate for amplification on the LightCycler 480 PCR system (Roche). The steps 
of the cycling program were 95 °C for 10 min followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. Data were 
analyzed using the online software provided by SABiosciences (http://pcrdataanalysis.sabiosciences.com/pcr/ 
arrayanalysis.php). 

Real-time RT-PCR. Total RNA in chicken lung was extracted as described above. Real-time RT-PCR was 
performed with QuantiNova Reverse Transcription kit (Qiagen, Alameda, CA) and QuantiNova SYBR Green 
PCR kit (Qiagen) using previously described primers that target the chicken TLR3, TLR4, and 28S 54 . All reactions 
were set up in duplicate on the iQ5 real-time PCR system (Bio-rad). The steps of the cycling steps were 95 °C for 
2 min followed by 40 cycles of 95 °C for 5 s and 60 °C for 10 s. Obtained cycle threshold (Ct) values were normal¬ 
ized to 28 S. The relative expression of TLR3 and TLR4 (fold change of naive control) was determined using the 

2 ~AAct me thod 55 . 

Statistical analysis. Data were analyzed by unpaired t tests or ANOVA followed by Dunnett’s multiple com¬ 
parison tests using GraphPad Prism (GraphPad Software, San Diego, CA). The p values smaller than 0.05 were 
considered significant. 
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